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Abstract 

A discrete complexified quaternion Fourier transform is introduced. Tliis is a generalization of the 
discrete quaternion Fourier transform to the case where either or both of the signal/image and the 
transform kernel are complex quaternion- valued. It is shown how to compute the transform using four 
standard complex Fourier transforms and the properties of the transform are briefly discussed. 

1 Introduction 

Discrete quaternion Fourier transforms have been described by several authors [2,5,10,12,14]. The pioneering 
work of Ell [3] predated all of these papers and was the inspiration for [12,14]. Quaternion, or at least vector- 
valued, signals arise in color image processing and in processing of signals captured by vector sensors, for 
example, vector geophones [9]. The main motivation for using quaternion algebra to handle such signals 
is the correspondence between vector samples and the vector parts of quaternions, which permits holistic 
processing of the vector samples and signals, making use of information about the direction and magnitude 
of the samples in sample space. 

Since the introduction of the quaternion Fourier transform there has been a concern that perhaps a 
quaternion-valued signal or image requires a transform based on a 'higher' algebra, just as the Fourier 
transformation of real-valued signals or images requires a complex transform. Indeed Biilow and Sommer [2] 
have shown the advantages of using a quaternion transform even for real (grey-level) images, in order to 
analyse symmetries in the image. Unlike Biilow and Sommer, we are interested in signals and images with 
vector samples (that is, three components at least, perhaps even three complex components). 

However, until now there have been no reported definitions of Fourier transforms for quaternion-valued 
signals or images based on an algebra of higher dimension than the real quaternions*. 

At various times the octonions have been suggested as the basis for a higher transform, but their algebra 
is non-associative and therefore presents formidable problems compared to the quaternion algebra which 
is of course non-commutative, but otherwise relatively straightforward to work with. Therefore we have 
considered the case of a transform based on complexified quaternions, that is quaternions with four complex 
components. These can be considered equivalently as complex numbers with quaternion real and imaginary 
parts. The properties of this and other quaternion algebras have been extensively studied in the abstract 
algebra community, but much less knowledge is available at the practical algebraic level needed for engineering 
use. The complexified quaternions were studied by Hamilton himself [8] and [7, Lecture VH, §669 (p. 664)], 
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but he called them biquatemions. We prefer to avoid this term as it was also used by Clifford with a different 
meaning. 

fn this paper we introduce the bare basics of complex quaternions necessary to the paper and we define 
a complex quaternion Fourier transform analogous to the quaternion Fourier transforms discussed in [14]^. 
We then show how the transform may be computed, fn common with real quaternion Fourier transforms, 
a direct (DFT or FFT) implementation is possible using a complex quaternion arithmetic package but it is 
also possible to factorize the transform into four complex transforms, and thus compute it using standard 
complex FFT code, thus exploiting the significant effort expended on performance enhancement of FFT code 
by other authors. 

The details of quaternions and of Fourier transforms in general are outside the scope of this paper. 
Ward's book [15] provides a good introduction to quaternions, and there are many texts in the area of 
signal processing that discuss discrete and fast Fourier transforms. Generalizing the quaternions to complex 
quaternions is straightforward - the four components of the quaternion take complex values rather than 
real values, but there are some important ramifications. Firstly, although the quaternions are a division 
algebra, the complex quaternions are not. Any real quaternion (other than zero) has a multiplicative inverse 
and a non-zero norm or modulus. In contrast, a complex quaternion (other than zero) may have a zero 
norm or modulus and therefore no multiplicative inverse. Obviously this means that care is needed when 
computing with complex quaternions, and if we are to define and compute Fourier transforms with them, 
we must establish whether vanishing norms will cause problems. A second issue is that a quaternion Fourier 
transform necessarily contains as its kernel a quaternion exponential function, and if we are to generalize 
this to complex quaternions, we need a complex quaternion exponential function. As has been shown in [14], 
quaternion exponentials can be defined using an arbitrary unit pure quaternion since these are roots of — 1. 
(It is fundamental to Fourier transforms that a root of —1 is needed, since the transform analyses a signal 
into sinusoidal components, and the famous formula of de Moivre: e'^ = cos^ + isin^ generalizes to any 
algebra in which a root of —1 can be defined. This follows from the series expansions of the exponential, 
cosine and sine functions.) If we are to define a complex quaternion Fourier transform, we require a complex 
quaternion root of —1. We cannot use a real quaternion root of —1, or the transform will simply be a 
quaternion Fourier transform even if it operates on a complex quaternion signal (the transform would be 
equivalent to two quaternion transforms applied separately to the real and imaginary quaternion parts of the 
signal). The solutions of the equation q"^ = —1 [13] are therefore crucial to defining the transforms presented 
in this paper. 

2 Definitions 

Before defining the transform we must define some complexified quaternion concepts. We have not given 
proofs of these: most can be found in [15] although with different notation. A complex quaternion has four 
complex components, for example: q = w + xi + yj + zk, where w, x, y, z are complex and i, j and k are the 
three quaternion operators. We denote the complex square root of —1 by a capital I to distinguish it from the 
quaternion i with which it must not be confused (this is a fundamental aspect of complexified quaternions). 
Note that, although the multiplication of (complex) quaternions is not commutative, I, and therefore all 
complex numbers (complex scalars), commute with i,j and fc, and therefore with all quaternions. In what 
follows, 3?(g) denotes the real part of q, that is: ^{q) = 3?(w;) + ^{x)i + ^{y)j + ^{z)k, and similarly 9(g) 
denotes the imaginary part of q. 

2.1 Inner product, norm and modulus 

The inner product of two complex quaternions q = Wq + Xqi + yqj + Zqk and p = Wp + Xpi + ypj + Zpk is 
(<?) p) = WqWp + XgXp + yqyp + ZqZp. The semi-norm^ of a quaternion \ \q\ \ = {q,q) = w"^ + x'^ + y'^ + z"^ , and the 
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Here, it is possible for = 0, even though g ^ as discussed in subsection 2.3. 
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modulus is the square root of the norm. All of these results are complex, in general. In the special case where 
= the semi- norm reduces to the usual quaternion norm. If the quaternion has zero real part {^{q) = 0), 
the semi- norm reduces to a real, but negative, norm (because = —1). The inner product provides us with 
the concept of orthogonality, which is essential to the construction of a (complex) orthonormal basis. Two 
complex quaternions are orthogonal if their inner product is zero, that is g _L p ■^=^ {l^P) = (Pi ?) = 0- 

2.2 Complex quaternion roots of —1 

Let /X be a complexified quaternion root of —1. It was shown in [13] that /lx must be a pure complex 
quaternion satisfying the following conditions: 

5R(m)T3(m), ||5R(/x)||-||3(/x)|| = 1 
For example, n = i+ j + k + {j — fe)I is a root of —1 and has a unit (real) norm, as can easily be verified. 

2.3 Complex quaternions with null modulus 

It is possible for a non-zero complex quaternion to have zero semi-norm. The conditions for the norm to 
vanish were discovered by Hamilton [7, Lecture VII, §672 (p. 669)]. Informally, the result can be shown as 
follows. Starting with the semi-norm of the quaternion in the form \\q\\ = w"^ + + + , expand one of 
the components in terms of its real and imaginary parts, for example w: 

w'^ = (SR(w) + 3(w)I)2 

= ?ft{wf - Q{wf + 2SR(w)9(w;)I 

We see that the real part of the semi-norm consists of the difference between the norms of the real and 
imaginary parts of q, and the imaginary part of the semi-norm is twice the inner product of the real and 
imaginary parts of q. For the semi-norm to vanish, its real and imaginary parts must separately be zero. 
This requires the real and imaginary parts of q to be orthogonal real quaternions of equal norm, that is: 
(K(9),9(g))=0andp(g)|l = p(g)||. 

2.4 Transform pair 

Given the definition of the transform 'axis'§ as any complexified quaternion root of —1, /x, we can now define 

the transform itself, which is almost trivial. In this paper we restrict our account to the one-dimensional 
case, but all our definitions and results generalize to the two-dimensional case without difficulty. 

F[u] = J2 (-2^/^^) 

ra=0 
N-1 

/N = ^ ^exp(27rM^)FM 

u=0 

where the signal f[n\ and its 'spectrum' F[u\ have samples. The placement of the minus sign in the 
forward transform is, of course, arbitrary. The transform pair defined in equation 1 must, of course, be 
evaluated in complex quaternion arithmetic, and the exponential is a complex quaternion exponential. All 
the necessary operations, and the transform itself, have been implemented in the Quaternion Toolbox for 
Matlab@ library [11]. It is conventional in Matlab® for the scale factor to be applied to the inverse 
transform and we have followed this convention in implementing the complex quaternion transforms^. 

§The axis of the transform is a direction in 3-space (in this paper complex 3-space) which defines the orientation of the sine 
component of the transform in the space of the vector part of the signal samples. 

'An alternative and widely used convention is to distribute the scale factor between the forward and inverse transforms. If 
the two scale factors are to be equal, each must be 
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An alternative transform is possible, by interchanging the positions of the exponential and the function 
f[n\. This transform is closely related to the first and easily implemented, although we omit the details here. 

Now we consider whether vanishing norms can occur in the transform. The norm of the exponential 
cannot vanish for any n or u because fi has unit modulus, and the cosine and sine parameters are real. 
However, the signal f[n] could contain samples with vanishing norms, unless there is some restriction on 
the sample values which prevents the conditions for a vanishing norm occurring. For example, if the sample 
values are real or imaginary quaternions they cannot have vanishing norms. In the more general case where 
the samples have complex quaternion values, the possibility of samples with vanishing norm cannot be 
eliminated. Such samples would not contribute to the 'spectrum' and therefore could not be reconstructed 
by the inverse transform. It is not known whether samples of the spectrum could have vanishing norms 
if computed from real (or imaginary) quaternion signals. If this did occur, it would again mean that the 
transform would not correctly invert. Clearly this point requires further study. 

3 Factorization into four complex transforms 

The definition of the transform in (1) leads directly to a discrete Fourier transform implementation, which 
is useful as a reference; or test implementation, but useless for practical computation except for small values 
of A''. A fast implementation requires either a custom coding of a classic FFT algorithm in complexified 
quaternion arithmetic, or a decomposition into complex Fourier transforms that can be implemented using 
existing complex functions or code. It is the latter approach that we describe here and it follows from a 
similar approach developed for the real quaternion case in 2000 [4]. 

The signal function f\n] in (1) may be expressed in terms of an orthonormal basis defined by the 
transform axis. Apart from being a complex orthonormal basis, this is exactly as in the real quaternion case. 
The complex orthonormal basis is defined by /j, and two other unit pure complex quaternions v and ^ such 
that n ± f ± ^ and /xu = ^ (and /xf^ = —1)11. The basis may also be represented by a 3 x 3 complex 
orthogonal {not Hermitian) matrix: 



where /x = fixi + l^yj + i^-zk etc. We show here only the decomposition of the forward transform, as the 
inverse is similar. We first write /[n] in terms of its four complex components: 



Now applying a change of basis, we express f[n] in terms of the new basis. The change of basis is very simply 
implemented by resolving the vector part of each sample of f[n] into the three complex directions defined 
by the new basis, using the inner product defined in subsection 2.1. Thus we have: 




f[n\ 



w[n] + x[n]i + y[n]j + z[n]k 



x'[n] = (/x, x[n]i + y[n]j + z[n]k) 
y'[n] = {if, x[ri\i + y[n]j + z[n]k) 
z'[n\ = (4, x[n\i + y[n\j + z[n\k) 



In the new basis, we have: 



f[n\ = w[n] + x[n]ii + y'[n]i> + z'[n]$^ 
= [w[n] + x'[n]ij] + [y'[n] + z'[n]ij] v 



Using this change of basis, we are able to separate the transform into the sum of two transforms: 




\x i-v means (/l/., i^) = 
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Wc now separate the terms on the right into real and imaginary parts and group the real components together 
and the imaginary components together to make four complex terms. 

""t^ [5R(w[n]) + sR(x'[n])/i] \ 

"'^T? [»(t/'[n]) + 3fJ(z'[n])/x]j^ 
■'^^[9(y'[n])+3(^'[n])/i] i^J 

All four of the transforms within this expression are now isomorphic to a complex Fourier transform. That is, 
we may replace fj, with I (the complex root of —1, not the quaternion root of —1, i) in order to compute the 
transform, and we will obtain the same numeric results. After computing the four complex transforms, all 
that remains is to re-assemble the parts of F[u] and invert the change of basis. The latter step is equivalent 
to multiplying out the factors of fi and v> appearing above, but it is more easily performed by a change of 
basis using the transpose of the original basis matrix used to change from the standard basis to the (/x, i/, ^) 
basis. 

Note that the factorization into four complex transforms reduces to a factorization into two complex 
transforms in the special cases where the function to be transformed has cither zero imaginary part or zero 
real part. If, additionally, fi is restricted to be real, the transform reduces to a previously published [4] 
quaternion Fourier transform. However, if the function to be transformed is either real or imaginary, the 
transform defined in this paper is not equivalent to a real quaternion transform, because the kernel is a 
complex quaternion exponential. Finally we note that if the function to be transformed has no vector part 
(that is, it is real or complex), and we choose for /x the degenerate value I (the complex root of —1) the 
transform reduces to the usual complex Fourier transform. 

The factorization described above has been used in the implementation of the transforms described in 
this paper in both one and two dimensions as part of the Quaternion Toolbox for Matlab@ library [11, vO.8 
onwards]. A direct DFT implementation of equation 1 is provided for verification of the fast transforms. The 
complex fast transforms are the usual Matlab® functions, which are implemented internally by the FFTW 
library [6]. 

4 Discussion 

The complex quaternion transform applied to a real quaternion signal exhibits quaternion conjugate sym- 
metry in the 'spectrum' in the same way as a complex transform applied to a real signal. The shift and 
similar properties follow those of the complex transform. It remains to be seen in what ways the complex 
quaternion transform may be useful, but one possible area is the development and fast implementation of 
linear vector filters based on convolution with complex quaternion coefficients. We leave a more detailed 
discussion of these issues for a later paper. 

5 Conclusions 

This paper has introduced for the first time a complex quaternion Fourier transform which may be used to 
compute the transform of a quaternion or complex quaternion signal or image. It may be easily computed 
using decomposition into four complex Fourier transforms. It appears that the transform exhibits symmetries 
that are absent in quaternion Fourier transforms, just as the complex Fourier transform shows symmetries 
for real signals that are not apparent with a real-valued transform such as the Hartley transform or the 
DCT. 



N-l 



F[u] = 



rt=0 



V + e- 



5 



References 



[1] E. J. Borowski and J. M. Borwein, editors. Collins Dictionary of Mathematics. HarperCollins, Glasgow, 
2nd edition, 2002. 

[2] T. Biilow and G. Sommcr. Hypcrcomplcx signals - a novel extension of the analytic signal to the 
multidimensional case. IEEE Trans. Signal Processing, 49(ll):2844-2852, November 2001. 

[3] T. A. Ell. Hypcrcomplcx Spectral Transformations. PhD thesis. University of Minnesota, 1992. 

[4] T. A. Ell and S. J. Sangwinc. Decomposition of 2D hypcrcomplcx Fourier transforms into pairs of com- 
plex Fourier transforms. In Moncef Gabbouj and Pauli Kuosmanen, editors, Proceedings of EUSIPCO 
2000, Tenth European Signal Processing Conference, volume II, pages 1061-1064, Tampere, Finland, 
5-8 September 2000. European Association for Signal Processing. 

[5] M Felsbcrg and G Sommcr. Optimized fast algorithms for the quaternionic Fourier transform. In 
F. Solina and A. Leonardis, editors. Computer Analysis of Images and Patterns, volume 1689 of Lecture 
Notes in Computer Science, pages 209-216. Springer, 1999. 8th International Conference, CAIP'99 
Ljubljana, Slovenia, September 1-3, 1999 Proceedings. 

[6] FFTW. [Online], 1997-9. Software library, written by Matteo Prigo and Steven G. Johnson at Mas- 
sachusetts Institute of Technology. Available: http://www.fftw.org. 

[7] W. R. Hamilton. Lectures on Quaternions. Hodges and Smith, Dublin, 1853. Available online at Cornell 
University Library: http : / /historical . library . Cornell . edu/math/. 

[8] W. R. Hamilton. On the geometrical interpretation of some results obtained by calculation with bi- 
quaternions. In H. Halberstam and R. E. Ingram, editors. The Mathematical Papers of Sir William 
Rowan Hamilton, volume III Algcibra, chapter 35, pages 424-5. Cambridge University Press, Cambridge, 
1967. First published in Proceedings of the Royal Irish Academy, 1853. 

[9] Nicolas Le Bihan and Jerome Mars. Singular value decomposition of quaternion matrices: a new tool 
for vector-sensor signal processing. Signal Processing, 84(7):1177-1199, July 2004. 

[10] Soo-Chang Pei, Jian-Jiun Ding, and Ja-Han Chang. Efficient implementation of quaternion Fourier 
transform, convolution, and correlation by 2-D complex EFT. IEEE Trans. Signal Processing, 
49(ll):2783-2797, November 2001. 

[11] Quaternion Toolbox for Matlab®. [Online]. Software library, written by Stephen Sangwine and Nicolas 
Lc Bihan. Available: http://qtfm.sourceforge.net/. 

[12] S. J. Sangwine. Fourier transforms of colour images using quaternion, or hypercomplex, numbers. 
Electronics Letters, 32(21):1979-1980, 10 October 1996. 

[13] S. J. Sangwine. Biquaternion (complexified quaternion) roots of-1, June 2005, arXiv:math.RA/0506190. 

[14] S. J. Sangwine and T. A. Ell. The discrete Fourier transform of a colour image. In J. M. Blackledge 
and M. J. Turner, editors. Image Processing II Mathematical Methods, Algorithms and Applications, 
pages 430- 441, Chichester, 2000. Horwood Publishing for Institute of Mathematics and its Applications. 
Proceedings Second IMA Conference on Image Processing, De Montfort University, Leicester, UK, 
September 1998. 

[15] J. P. Ward. Quaternions and Cayley Numbers: Algebra and Applications, volume 403 of Mathematics 
and Its Applications. Kluwer, Dordrecht, 1997. 



6 



